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■ Abstract. In this paper I summarize the science motivations, as weU as a few 

' mid-infrared spectroscopic methods used to identify the principal mechanisms of 

^SJ , energy production in dust enshrouded galactic nuclei. The development of the 

various techniques is briefly discussed. Emphasis is given to the use of the data 
which are becoming available with the infrared spectrograph (IRS^) on Spitzer, 
as well as the results which have been obtained by IRS over the past two years. 
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Introduction 



^ ■ One of the open issues in extragalactic astronomy is how to quantify the phys- 

\ ical mechanisms that contribute to the energy production in galactic nuclei. 

■ The most luminous galaxies often display strong evidence of massive star forma- 

tion taking place in their nucleus. This is typically due to the large quantities of 
' atomic and molecular gas driven in their center as a result of instabilities on their 

. disk such as bars, which often form during tidal interactions. It is also widely ac- 

' cepted that m ost galactic nuclei likely harbor a super-massive black hole (SMBH) 

r~| . (see review of iFerrarese fc Fordll2nn,^ b Due to the high densities found in the 

Qh' nuclei the SMBH may accrete matter at a variable rate, even though this ac - 

cretion does not always result in emission of radiation ( Naravan et ahlflOQ^T ). 



I which would characterize the galaxy as harboring an active nucleus (AGN). A 

^ . prime example of such "non-AGN" is our own Galactic Center, which does con- 

tain a ^lO^M fT, SMBH and accretes material, but does not display any visible 
AGN activity ()Genzel fc Towneslll987l : iGhez et aLlligQ^ ). However, in the cases 



where AGN activity is present, hard electromagnetic radiation, mostly in form 
of UV and X-rays, is emitted from the nucleus. This is typically observed either 
\ directly, in the form of X-rays and radio emission which are less affected by 

obscuration, and/or indirectly as high energy photons are absorbed by the gas 
and dust and subsequently re-emitted at longer wavelengths. So the original 
question naturally translates to what fraction of the bolometric energy observed 
from a galaxy originates from an accretion disk (AGN) compared to usual star 
formation activity. This issue is of particular interest in cases where most of 
the energy of a galaxy originates from the nuclear regions rather than the outer 
regions / disk. 

Historically, our knowledge on the physics of the nuclear activity of galaxies 
is mostly based on optical and near-IR spectroscopy, obtained in major ground 
based telescope facilities. A wealth of diagnostic methods to classify and quantify 
the AGN activity have been developed over the years (ie lVeilleux &: Osterbrockl 



^The IRS was a collaborative venture between Cornell University and Ball Aerospace Corpora- 
tion funded by NASA through the Jet Propulsion Laboratory and the Ames Research Center. 
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119871 : lArmus et al.lll989l : lKennicutdlT993 : lKewlev et al.ll2nnih . A presentation of 
these methods is well beyond the scope of this review but it is only fair to state 
that they have been extremely successful in providing insight into the properties 
of galactic nuclei, in particular in cases where a direct line of sight to the nucleus 
is available. Polarization measurements have also been very useful in identifying 
evidence of AG N activity, in case s where the nucleus is indirectly probed via 
scattered light ()Heisler et al.lll997^ ■ 



The effect of dust absorption though, which does vary substantially as a 
function of wavelength, is one of the principal reasons which complicate the 
interpretation of the optical data. It is useful to remind ourselves that a lOeV 
optical photon can penetrate ~0.5mag of dust, while a IkeV X-ray photon can 
pass through a hydrogen column of ~ lO^^cm"^. If we were to consider the case 
of our Galactic Center where the Ay ~30mag, only ~1 in lO^'^ optical photons 
emitted in the nucleus of the Milky Way can reach our Sun. However, if we were 
to observe in K-band (2.2/im) the A2,2/im ~2.5mag, and one can detect 10% of 
the near-IR photons emitted from the source. 

Based on the arguments above. X-rays are clearly the best diagnostic of an 
AGN especially for higher redshift systems, as the slope of the X-ray spectrum 
and th e sensitivity of the detec tors result in relatively flat sensitivity out to 
z~4-5 ([Brandt Sz Hasingedl20o3 ). However, the energy emitted in X-rays from 
the majority of luminous extragalactic sources is ~ lO^'^ergs"^, several orders 
of magnitude bellow their bolometric luminosity. This, in addition to the fact 
that X-ray observations need to be performed from space, and most extragalac- 
tic sources are intrinsically faint in X-rays, makes estimates of the contribution 
of an AGN to their bolometric luminosity rather challenging (see iMushotzkvl 
iool). Radio waves are an appealing alternative, in particular given the latest 
improvements in VLBA ((Lonsdale et al. 1 1199.1 l2nn.il ^. Even in this case though, 
there are difficulties in detecting higher redshift systems, as well as in disentan- 
gling issue s related to no n-thermal emission and self-absorption from embedded 
nuclei f see ICondonlll99 ^ . 

The use of the infrared part of the spectrum 3 < A < lOOO^m in order 
to develop AGN diagnostics is an avenue which has been explored for the past 
twenty years for two main reasons: a) most luminous galaxies, QSOs being a 
"bright" exception, emit a sizeable fraction their energy (>30% in the infrared 
and b) the infrared is considerably less affected by dust absorption than the 
optical. 

The all sky survey by IRAS provided the first opportunity to classify a large 
numbe r of galaxies based on t he shape of their global spectral energy distribution 
fSEDl. lde Griip et al.l (|l987j 'l showed that a ratio of S60/S25>0.262 could result 
in a 70% success rate in identifying previously unknown Seyferts. More recent 
developments on broadband mid-IR diagnostics are presented in another review 
on this volume. 

2. Mid-IR Emission Line Diagnostics 

The advent of mid-IR spectroscopy with the Infrared Space Observ atory (ISO) 
opened a new window in the study of active and starburst galaxies fsee lVerma et al.1 



^Where S25 and S60 are the 25 and to 60 fim IRAS flux densities. 
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Figure 1. Mid-IR excitation diagrams for [NeV]A14.3/im (top) and 
[OIV]A25.9yLtm (bottom) as a function of the 6.2/^m PAH EW, for starbursts 
(green circles - Brandl et al. 2006), AGN (red squares, Weedman et al. 2006) 
and ULIRGs (in blue diamonds), adapted from Armus et al. (2006). The 
vertical (red) and horizontal (green) lines indicate the fraction of AGN and 
starburst contribution. 

l2fln.'Tl . for a review). Emission lines from relatively high excitation ions which 
are not produced by O stars can be used as a direct proble of AGN activity. 
In the infrared such lines are the [NeV]A14.3/xm and [NeV]A24.3;um with an 
ionization potential of 97.1eV, and the [OIV]^A25.9;um with 54.9eV. The neon 
hues were detected by ISO in Seyferts, though their detection in ULIRGs was 
challenging du e to the limitations in the sensitivity of the ISO spectrographs 
( Genzel et al.lll99Hl : ISturm et 31.1 12000') . Another indicator is the strength of the 
PAH emission bands which are suppre ssed in AGN sinc e the hard radiation field 
photo-dissociates the PAH molecules ()Lutz et al.lll998t ). 



■^The [OIV] line was easier to detect by ISO, but its predictive power is not as strong since it 
can also be produced in areas of strong galactic outflows and shocks. 



4 



Charmandaris 



The improved sensitivity of IRS on Spitzer ( Houck et al .11200^ : 1 Werner et alJ 



a factor of ~100 compared to ISO, enables us to extend these diagnostics. 
More detections based on the [NeV] hnes were secured for a large sample of 
galaxies and an in depth study of PAH emission has become possible (see Fig. 1). 
Moreover, techniques employing correlations based on the strength of other lines 
such as [SIV]A10.5/im an d [SIII]A18.7/33.5A^m or pricipal cornponent analysis are 
also being explored (see fWeedman et al.ll2005l : lArmus et al. I l2nn(il : iBrandT et al.l 
\20()(i : iDale et al.ll200g : lBuchanan et al.ll2006l . for details). More specifically the 
use of the 6.2;^m PAH is now readily used instead of the 7.7fim PAH which is 
more affected by uncertainties in the extinction due to the absorption by the 
9.7/um silicate feature. 

In addition to the mid-IR high-ionization lines and the 6~18/im PAH emis- 
sion features, the strength of the 3.3/im PAH has also been proven a useful AGN 
tracer. This feature, which can be probed from the ground for relatively bright 
sources, h as an EW of ^ 0.1/im in starburst systems but is also supressed in 
AGN fsee llmanishlll2nn(ir) . 

3. Mid-IR Continuum Diagnostics 

Even though high ionization line emission is the ideal probe of an active nu- 
cleus, despite the IRS sensitivity, the lines are often difficult to detect in highly 
extinct and faint/distant sources. The presence of an AGN though can also 
be inferred by detecting the thermal emission from dust surrounding the pu- 
tative torus heated by the accretion disk to nearly sublimation temperatures 
(~1000K) and radiating in near equilibrium. This emission had been detected 
in Seyfert galaxies and is now evident in IRS spectra of distant sources (see 
lAlonso-Herrero et al.ll200^ . and Fig. 2). 




Figure 2. The SED of the nucleus of NGC1068, our nearest Sy-2 at 
D=14.4Mpc (left), adapted from Alonso-Herrero et al. (2003) compared to 
the IRS spectrum of FSC15307+3252, a quasar at z=0.92, D=5.8Gpc (see 
Teplitz et al. 2006). Note the similarity in the thermal emission from the 
dust near the AGN torus heated to near sublimation temperatures. 

Unlike extreme starbursts or HII regions which also display a rising slope 
at A > Wfim, due to the heating of the grains to T~300K by the embedded 
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0/B stars, the pres e nce o f this hot continuum at A ~ 5/im is only seen in 
AGN. iLaurent et alJ (|2000l ) proposed an AGN diagnostic method - now called 
the "Laurent Diagram" - which takes advantage of this difference, as well as the 
destruction of PAHs in both AGN and extreme starbursts. Using three template 
mid-IR spectra for a "pure" AGN, an HII region, and a photodissociation region 
(PDR) they defined an AGN dominated locus in a two parameter phase-space 
(see Fig. 3). The proximity of a galaxy to one of the three corners, i.e. the AGN 
corner, would suggest the extend by which the mid-IR spectrum displays AGN 
characteristics. 




Figure 3. Application of the "Laurent Diagram" using IRS spectra. On 
the left we display the three template spectra: M17 (scaled by 10~^) as the 
Hll/extreme starburst, 3C273 as the pure AGN template, and NGC7023 as 
the pure PDR template. On the right we show the location of a number of 
starbursts, AGNs and ULIRGs on the diagram. The dotted and solid lines 
indicate the 50 and 75% contribution of the corresponding template to the 
integrated spectrum of a source (see also Armus et al. 2006) 

A first application of this m ethod to Spitzer/IRS spectroscopy has been 
presented in I Armus et al. In Figure 3 we display the principle of the 

method with the three template spectra used, as well the classification of a 
number of galaxies based on it. Even though we suggest the reader to review 
the original paper for a detailed discussion of the method, a few points are worth 
mentioning here. 

• The predictive power of the method depends strongly on the selection of 
the template spectra used as the three cornerstones. In the original paper 
of lLaurent et al.. (2000, 1 Centaurus A was selected as the AGN template, 
while now the wealth of high quality IRS spectra enabled us to fine tune 
the selection. 

• The actual coordinates of the three templates depend on the corresponding 
PAH to 5. 5/im continuum and 14.5 to 5. 5/im continuum ratios. Since the 
continuum fluxes are calculated by integrating under the spectrum, the 
wavelength range used for the integration directly affects the value of the 
ordinate for all points in the diagram. It is not obvious that a flat ufi, 
spectrum will result in an ordinate equal to 1. In addition, the exact 
selection of the limits and underlying continuum in the calculation of the 
6.2/im PAH, also affects the value of the abscissa. 
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• None of the templates have been corrected for extinction, so the intrinsic 
extinction in the observed integrated spectrum of a source may result in 
placing it outside the parameter space of the diagram. Note that the 
nucleus of a galaxy, which may contain a strong AGN and/or HII region 
component, is often more enshrouded than the PDR regions in its disk. 
As a result of this variable extinction, there is no one-to-one translation 
of an extinction vector for all points placed in the Laurent diagram. 

• The method provides quantitative means for estimating the AGN contri- 
bution to the "integrated mid-IR spectrum" of a galaxy. Even though the 
~ 25-35/um flux traces luminosity it has not yet been demonstrated that 
the this method can be used to estimate the contribution of an AGN to 
the infrared luminosity (Lir) of a system. 

The "Laurent Diagram" appears to be the most powerful method for AGN 
classification in the mid-IR for galaxies where no detections or strong limits on 
high ionization lines are available. More analysis based on Spitzer/IRS spectra 
will be need to extend its predictive power to the whole IR. Ongoing work 
indicates that the method can also be applied in Spitzer/IRS low resolution 
spectroscopic studies of extended star forming regions where variations in the 
slope of the mid-IR spectrum and the strength of PAHs can be used to quantify 
the intensity of star formation of the regions (Leboutellier priv. comm.). 

4. Theoretical Modeling 

Another approach to quantitatively assess the contribution of the various com- 
ponents in the observed IR emission from a source is to fit the IR spectrum using 
theoretical models. Ideally one could use a full 3D radiative transfer calcula- 
tion to model the dust properties and match the observed spectrophotometry. 
A large set of t heoretical SEDs are now available for comparison with th e ob- 
servations (i.e. lElitzur Sz Shlosmanl l2006l : ISiebenmorgen Sz Krueeell l2006l . and 
references therein.). However, limited knowledge of the details in the spatial 
distribution and geometry of the sources, in particular for systems such as star- 
burst and ULIRGs, which are inherently disturbed, often make this approach 
rather challenging. 

Recently, a new fitting approach relying on the Spitzer/IRS 5-38^m spec- 
tra, using constrains from near- and far-IR observations has been developed 
( Marshah et~ani2nnfir ). The method assumes that the SED of a galaxy can be de- 
composed into dust components at different characteristic temperatures, source 
emission components, embedded photospheric emission from starburst cores or 
an active nucleus as well as PAH feature emission. The emission from each 
component is calculated using a realistic dust model consisting of a distribution 
of thermally emitting carbonaceous and silicate grains. The model accounts for 
stochastic emission from very small grains by fitting model PAH templates to 
the spectra. This method has been applied in a variety of Spitzer/IRS data and 
most sources are well fit by a maximum of four components referred as cold 
(~30K), tepid (~100K), warm (~200K), and hot (>300K). The results of the 
fits are used to calculate the relative bolometric contributions from the different 
components, providing a method to compare the infrared properties of star- 
burst galaxies and quasars. Sample results of the fitting method on NGC7714 - 
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a starburst galaxy - a ULIRG such as NGC6240 and the quasar PG0804+761 
are presented in Figure 4. 




Figure 4. Model fits to tlic IRS spectrum and SED of NGC7714, NGC6240, 
and PG0804+761. The contribution of the various temperature components 
is indicated by the dotted and dashed lines. Note how well the model fit trace, 
PAH emission features, as well as the silicate absorption and - in the case of 
PG0804+761 - emission features. The presence of an AGN in PG0804+761 is 
also revealed by the hot dust component peaking at ~ 4^m. For more details 
on the luminosity contribution of each component see Marshall et al. (2006) 



5. Conclusions 

It has become evident from the numerous contributions presented during this 
meeting that the high quality of Spitzer/IRS spectroscopy is opening new hori- 
zons in the use of the infrared as a tracer of the properties of nearby and high- 
redshift galaxies. With more than three years of mission to go, we have just 
glimpsed on the possibilities that lie ahead. 



8 



Charmandaris 



Acknowledgments. I would like to thank Jim Houck for giving me the op- 
portunity to spend 5.5 formative years in Ithaca and participate in the Spitzer/ 
IRS "venture". I would also like to acknowledge all members of the IRS in- 
strument team at Cornell and Caltech for many enlightening and stimulating 
discussions on spectroscopy and the infrared, as well as the organizers for the in- 
vitation to this exciting conference. I also appreciate the help of Lee Armus and 
Jason Marshall who provided data and plots for this paper prior to publication. 



References 

Alonso-Herrero, A. et al. 2003, AJ, 126, 81 

Armus, L., Heckman, T. M., & Miley, G. K. 1989, ApJ, 347, 727 

Armus, L. et al. 2006, ApJ, (submitted) 

Brandl et al. 2006, ApJ, (submitted) 

Brandt, W. N., & Hasinger, G. 2005, ARA&A, 43, 827 

Buchanan, C. L., et al. 2006, AJ, 132, 401 

Condon, J. J. 1992, ARA&A, 3 0, 575 

Dale, D. 2006, ApJ, (in press - astro-ph/0604007) 

de Grijp, M. H. K., Lub, J., & Miley, G. K. 1987, A&AS, 70, 95 _ 

Elitzur, M., & Shlosman, I. 2006, ApJ, (submitted - astro-ph/0605686) 

Ferrarese, L., & Ford, H. 2005, Space Science Reviews, 116, 523 

Genzel, R., & Townes, C. H. 1987, ARA&A, 25, 377 

Genzel, R., et al. 1998, ApJ, 498, 579 

Ghez, A. M., Klein, B. L., Morris, M., & Becklin, E. E. 1998, ApJ, 509, 678 

Heisler, C. A., Lumsden, S. L., & Bailey J. A. 1997, Nat, 385, 700 

Houck, J. R., et al., 2004, ApJS, 154, 18 

Imanishi, M. 2006, AJ,131, 2406 

Kennicutt, R. C., Jr. 1998, ARA&A, 36, 189 

Kewley L. J., et al. 2001, ApJ, 556, 121 

Laurent et al. 2000, A&A, 359, 887 

Lonsdale, C. J., Smith, H. J., & Lonsdale, C. J. 1993, ApJ, 405, L9 

Lonsdale, C. J., Lonsdale, C. J., Smith, H. E., & Diamond, P. J. 2003, ApJ, 592, 804 

Lutz, D., et al. 1998, ApJ, 505, L103 

Marshall, J. A., et al. 2006, ApJS, (submitted) 

Mushotzky, R. "How are AGN Found?" in in " Supermassive Black Holes in the 
Distant Universe" (2004), ed. A. J. Barger, Kluwer Academic Publishers, 



astro-ph /0505144 

Narayan, R., et al. 1998, ApJ, 492, 554 

Siebenmorgen, R., & Kruegel, E. 2006, A&A, (submitted - |astro-ph/0606444 | 

Sturm, E., et al. 2000, A&A, 358, 481 

Teplitz, H.L, et al. 2006, ApJ, 638, 1 

Veilleux, S., & Osterbrock, D. E. 1987, ApJS, 63, 295 

Verma, A., Charmandaris, V., Klaas, U., Lutz, D., & Haas, M. 2005, Space Science 

Reviews, 119, 355 
Weedman, D. W., et al. 2005, ApJ, 633, 706 
Werner, M., et al., 2004, ApJS, 154, 1 



